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ABSTRACT 

Based on fractional Brownian motion (fBm) simulations of 3D gas density and 
velocity fields, we present a study of the statistical properties of spectro-imagery 
observations (channel maps, integrated emission, and line centroid velocity) in 
the case of an optically thin medium at various temperatures. The power spectral 
index of the integrated emission is identified with that of the 3D density field 
(7„) provided the medium's depth is at least of the order of the largest transverse 
scale in the image, and the power spectrum of the centroid velocity map is found 
to have the same index jc as that of the velocity field (7?,). Further tests with 
non-fBm density and velocity fields show that this last result holds, and is not 
modified either by the effects of density-velocity correlations. A comparison is 
made with the theoretical predictions of Lazarian & Pogosyan (2000). 



1. Introduction 



The early works on interstellar turbulence, summarized by Chandrasekhar (1949), pointed 
out that the interstellar medium (ISM) exhibits Reynolds numbers so large that it is very 
likely to be turbulent. For incompressible turbulent flows the Kolmogorov (1941) theory 
predicts a dissipation-free energy cascade between large scales where turbulent energy is 
injected, down to small scales where it is dissipated and heats the gas. In this inertial range 
the theory predicts a three-dimensional power law energy spectrum {E{k) (x fc~^/^), which 
has been observed extensively in terrestrial turbulence experiments (see Grant et al. (1962) 
for an example). 
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Since the work of Chandrasekhar (1949), several studies have shown that the interstellar 
medium density and velocity structures projected onto the plane of the sky are self-similar 
and well described by power laws. For instance, a number of works were done on Ha 
centroid velocity fields in Hll regions using auto-correlation and structure functions (see 
Miville-Deschenes et al. (1995) and references therein). The self-similar structure of Hi has 
also been studied, almost exclusively using the power spectrum of the 21 cm emission in 
galactic (Baker 1973; Crovisier & Dickey 1983; Green 1993; Deshpande et al. 2000; Dickey 
et al. 2001) and extra-galactic regions (Spicker & Feitzinger 1988; Elmegreen et al. 2001; 
Stanimirovic & Lazarian 2001). Various studies aimed at describing the fractal properties 
of dust infrared emission were also done using power spectra (Gautier et al. 1992; Wright 
1998) and the arca-pcrimctcr relation (Bazell & Desert 1988; Dickman et al. 1990; Vogelaar 
& Wakkcr 1994). But it is probably on molecular clouds (CO emission) that most of the 
work has been done using several statistical tools: the size-linewidth relation (Larson 1981; 
Myers 1983; Falgarone & Phillips 1990), the auto-correlation function (Kleiner & Dickman 
1984, 1985, 1987; Perault et al. 1986; Kitamura et al. 1993; Miesch & Bally 1994; Larosa 
et al. 1999), wavelet decomposition (Gill & Henriksen 1990), the area-perimeter relation 
(Falgarone et al. 1991), principal component analysis (Heyer & Schloerb 1997; Brunt & 
Heyer 2002b) and the A- variance method (Stutzki et al. 1998; Bensch et al. 2001). 

All those studies suggest that turbulence, for which such self-similar properties are ex- 
pected, may play a key role in the structuration of the interstellar medium, both in density 
and velocity. Furthermore, by providing non-thermal support against self-gravitational col- 
lapse, turbulence also regulates star formation (Klessen et al. 2000; Heitsch et al. 2001). It 
might also have a strong and very non-linear impact on dust evolution by accelerating the 
coagulation or fragmentation of grains (Falgarone & Pugct 1995; Miville-Deschenes et al. 
2002), and on the ISM chemical evolution as intermittent energy dissipation could locally 
trigger endothermic reactions (Joulain et al. 1998). However, following Scalo (1987), in- 
terstellar turbulence may be very different from Kolmogorov's prediction. The interstellar 
medium is compressible and magnetized, and several processes inject energy at very different 
scales (from galactic differential rotation down to stellar winds). Therefore, the description 
of interstellar turbulence in various environments is mandatory to understand its role on the 
structure, kinematics, energy transfer, thermodynamics and chemistry of the gas at various 
scales. 

One important issue here is to relate the statistical properties of the observed quantities 
(integrated emission, channel maps and centroid velocity) to the three-dimensional properties 
of the density and velocity fields, which actually describe interstellar turbulence. This calls 
for reliable methods to extract the relevant information from spectro-imagery observations, 
called position-position-velocity (PPV) data cubes, and disentangle the density, velocity, and 
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temperature contributions to the observed space-velocity structures. There is also a need 
to understand the effects of 3D-2D projection and of line opacity on the observed statistical 
properties. Lately there has been a number of works based on PPV data cubes to deduce the 
statistical properties of the underlying three-dimensional density and velocity fields of Hi 
observations (Stanimirovic & Lazarian 2001; Elmegreen ct al. 2001; Dickey et al. 2001). Most 
of these works are based on the theoretical work of Lazarian & Pogosyan (2000) (hereafter 
LPOO), who proposed a method to deduce the 3D density and velocity power spectral indexes, 
respectively noted 7„ and 7^ thereafter, from the PPV cubes. 

In this paper, we study the spectral properties of simulated PPV cubes using fractional 
Brownian motion (fBm) density and velocity fields, with the twofold purpose of comparing 
them with the predictions of LPOO, and proposing another method, based on the centroid 
velocity map, to determine the spectral index of the three-dimensional velocity field. The 
reason for using fBms, despite their lack of physical reality, lies with their well-behaved 
statistical properties and the ease with which they can be generated. With single-index 
power law power spectra and Gaussian fiuctuations, fractional Brownian motions represent 
useful and fiexible test cases for structure statistics retrieval methods. Such random fields 
were used for instance by Stutzki et al. (1998) to describe molecular cloud images using 
A-variance analysis, and by Brunt & Heyer (2002a) to retrieve the statistics of turbulent 
interstellar velocity fields through the use of Principal Component Analysis (PGA). It should 
however be stressed that these fields do not arise from physical simulations. In particular, 
they do not comply with hydrodynamics nor gravity, and are used solely for their statistical 
simplicity. 

This kind of study based on fBms usefully complements those conducted on magneto- 
hydrodynamical (MHD) simulations (see for instance Ostriker et al. (2001)) because the lack 
of physics is somewhat compensated by the broad range of scales accessible to the analysis 
in terms of power laws. In numerical simulations of turbulence the inertial range is severely 
limited by the available resolution and the large number of scales affected by injection and 
dissipation. 

The paper is organised as follows: In section 2, we briefiy present fBms and their proper- 
ties, in the context of designing simple model distributions for our three-dimensional density 
and velocity fields. The computation of PPV data cubes from these fields is detailed in sec- 
tion 3, then an analysis of their spectral properties is done in section 4. Section 5 presents 
our concluding remarks. An appendix contains some details on analytical computations. 



-4- 



2. Density and velocity fields 

2.1. Definition and construction of fBms 

A fractional Brownian motion (fBm) is a Gaussian random field F which can be defined, 
as done in Voss (1985), in any Euchdean dimension N by the relation 

< [F(r2) - F(ri)]2 > oc ||r2-ri|p^ (1) 

where H ^ called the Hurst exponent, is a real number in [0, 1], and the brackets stand for a 
spatial average over positions ri and r^- The value H = 0.5 gives the usual Brownian motion, 
for which it is well known that the mean squared increment < AF^ >=< [F{r2) — F(ri)]^ > 
scales as the separation ||Ar|| = ||r2 — ri||. This relation states that the fluctuations of the 
fleld are isotropic, and that their amplitude as a function of scale is a power law. It easily 
translates into a similar relation for the power spectrum, which is the Fourier transform of 
the autocorrelation function Af{p) —< F{r)F{r + p) >: 

Ppik) = ^^(k) oc A;^ (2) 

where isotropy is expressed by the fact that Pp depends only on the wavenumber k, which 
is the length of the A'"-dimensional wave vector k. The power spectral index 7 is related to 
the Hurst exponent and the dimension hy ^ — —2H — N . 

Several methods are available to numerically generate fBms (for instance as a collection 
of clumps with a power law mass spectrum, see Stutzki et al. (1998)), but the easiest is based 
on Eq. 2 and the fact that Pp is proportional to the squared amphtude of the Fourier 
transform of F. Hence we first compute an isotropic amplitude following a power law: 

A{k) = ^F/^ (3) 

where Aq is a normalization factor. Then the phase 0(k) is constructed using a random 
number generator algorithm based on a uniform distribution between — tt and tt, and to 
make sure the image in direct space is real, we enforce the constraint 0(— k) = — 0(k). From 
the amplitude and the phase, the real and imaginary parts of the Fourier transform are 
computed: 

Re[F(k)] = A{k) cos[0(k)] (4) 
Im[F(k)] = A{k) sin[(/.(k)] (5) 
and the fBm F(r) is simply obtained by taking the inverse Fourier transform of -F(k). 

The method described here is general for any dimension, and the power spectra of the 
fields generated this way obey exact power laws^, as can be seen in Fig. 1, which shows a two- 



The only random elements are the phases. 
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dimensional (2D) 257x257 fBm image, along with its power spectrum^. In this figure we used 
a scatter plot (all points of the 2D amplitude are shown) in order to assess the field's isotropy. 

But in other figures of the paper, where we compare power spectra, azimuthal averages are 
used to make the plots easier to read. It should also be noted that this Fourier-based method 
of generating fBms leads to periodic distributions and therefore our computations, unlike 
real observations, will not suffer from non-periodic boundary conditions^. 

2.2. 3D density and velocity fields 

The purpose of this paper is to study the effects of projection from a three-dimensional 
dynamical medium onto a spcctro-imagery observation. To conduct such a study, the first 
step is to generate model interstellar 3D density and velocity fields and from these compute 
simulated PPV data cubes. 

We therefore generated three-dimensional'^ fBms with spectral indexes of —3.0, —3.5, 
—4.0, and —4.5^. For each spectral index, two fBms were computed, one for the density field 
n{x,y,z) and one for the velocity field v{x,y,z). The computations were also performed 
with a uniform density field, for which the power spectral index is formally equal to —oo. 
Here we just consider the line-of-sight (longitudinal) velocity component, as it is the only 
one available to observers. Like Kitamura et al. (1993) we make the assumption that the 3D 
velocity field v{x, y, z) is isotropic and homogeneous and that the longitudinal and lateral 
velocity components have the same power spectra as v{x, y, z). 

For the sake of reahsm, some sort of correlation should exist between the 3D density and 
velocity fields. However, for simplicity, the main part of this paper deals with uncorrelated 
cubes, the effects of a certain amount of correlations between n and v being reported on in 
section 4.2. 

A typical line of sight extracted from one set of 3D density and velocity fields is shown 



^The odd size is simply a matter of convenience when building the fields. The Fourier transform algorithm 
is a bit slower, but speed was not critical in our simulations. 

^The reason why this is important is that numerical Fourier transform algorithms such as the FFT treat 
images as if they were periodic. If they are not, artificial jumps at the edges result in high frequency 
components which alter the power spectra. On the other hand, apodization of the images to prevent this 
may in turn affect large scales. 

*The dimensions are 129 by 129 by 129 pixels. 

^These correspond to Hurst exponents respectively equal to 0, 0.25, 0.5, and 0.75. 
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in Fig. 2. The density and velocity fields used in this example have power spectral indexes 
of 7^ =-3.0 and 7„ =-4.5 respectively. The former displays more small-scale structures, as 
expected for a shallower spectral index (7^ < 7„)^. By construction, the mean values of 
the original cubes are set by the normalization Aq and their variances are determined by 
the total number of pixels. To obtain physically plausible numbers, the following procedure 
was applied: The minimum values of the density cubes were subtracted from them in order 
to obtain only positive values, and a multiplicative factor was applied to obtain an average 
density of 90 cm~^. For the velocity cubes, we first subtracted the mean, then scaled the 
values to obtain a velocity dispersion of 3 km s~^. An important feature to point out is that, 
since typically only 0.1% of the points in a Gaussian distribution are more than 3a away 
from the mean, our density fields present quite low contrasts, with cr/n < 0.3. This is an 
unavoidable limitation of these model fields. 



3. Simulation of a PPV data cube 
3.1. Definition 

To obtain a PPV data cube from the simulated 3D density and velocity fields, we 
compute, at each position {x,y), the spectrum that would be observed in this direction. 
This is done under the simplifying assumption that the medium is observed in an optically 
thin line, such as, in some instances, the 21 cm line of neutral atomic hydrogen. In this 
case, the spectrum at a given position {x, y) is computed by simply adding the contributions 
of all cells along the fine of sight, the emission from a single gas cell being a Gaussian fine 
centered at velocity v{x,y,z) and of velocity-integrated area proportional to n{x,y, z)5z ^. 
Eventually, the PPV data cube Nu{x, y, u), which is the column density of Hi along the line 
of sight {x,y) at the velocity u within Su, can be expressed by the following equation^: 

/ ^ n{x,y,z)Sz ( [u - v{x,y,z)f\ 
N^{x,y,u) = l^ r—, rexp — ■ 6 



^This criterion holds when comparing fields of the same Euclidean dimension, but the relevant number in 
direct space when assessing the amount of small-scale structures is the Hurst exponent, as will be discussed 
later. Here 7„ < 7t, is equivalent to > 

''Here 5z is the depth, in cm, of a cell in the 3D density cube. The total depth of the cloud was set to 1 
pc and the size of the fBm is 129 x 129 x 129. In this case, = 2.414 x 10^^ cm 

^Here we consider the particular case where the spectral resolution of the observation is equal to the 
channel width and where the spectral transmission is a step function of width 5u. 
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The dispersion of the Gaussian is given by: 

(^(x,y,z) = 



dv{x,y,z) 



dz 



5z 



m 



(7) 



where ks is the Boltzmann constant, T is the gas kinetic temperature and m is the mass 
of the emitting species (in our case Hi). This broadening includes not only the thermal 
velocity dispersion ath = y^ksT/m, but also a small-scale velocity gradient which roughly 
models the subpixel structure of the velocity field (Brunt & Heyer 2002a). Also, to make 
sure that we lose a negligible fraction of the signal, the velocity limits in the final PPV cube 
are given by Vmm — 2crt/i and I'max + Scxt/i, where Wmin and Wmax are the minimum and maximum 
values in the 3D velocity field. The number of channels is then obtained by dividing this 
velocity range by the channel width Su, which was set to 0.25 km s"^. As a result, all of our 
PPV cubes do not have the same extent in the third dimension, since it depends on the 3D 
velocity cube and on the temperature. 



3.2. Synthesized spectra 

Examples of the fine synthesis procedure described above are shown in Fig. 3 for the 
same line of sight that was used in Fig. 2. These spectra were computed using equations 6 
for T —1 K, 10 K, 100 K, and 500 K. One can see the expected smearing effect of the 
temperature, which smooths the spectra as T increases, similar to the effect of decreasing the 
spectral resolution. As the temperature increases, the spectrum becomes smoother and closer 
to Gaussianity, indicating that the thermal broadening is close to the original dispersion in 
the input 3D velocity cube. For reference, the thermal dispersions for all temperatures are 
0.09 km s"^ at T =1 K, 0.29 km s"^ at T =10 K, 0.91 km s"^ at T =100 K and 2.0 km s"^ 
at T = 500 K. 

The spectrum computed at T = 100 K, a value representative of Hi gas, is reminiscent 
of what is observed at 21 cm in high latitude clouds. At artificially low temperature^ (T = 
1 K) many components appear. These components do not correspond to any physical "cloud" 
structure, as they simply represent the fractal properties of velocity and density fiuctuations 
along the line of sight. At such a low kinetic temperature, the line is far from the optically 
thin regime and the effect of opacity on the observed spectrum would be important. These 
spectra are shown only to illustrate the effect of kinetic temperature on our ability to resolve 



^In Hi, temperatures seem to vary from around 80 K in the cold neutral medium (CNM) to 5000 K - 
8000 K in the warm neutral medium (WNM). 
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kinematical structures. The effective spectral resolution of such PPV observations therefore 
includes the thermal broadening. Considering, as a simplification, a Gaussian-shaped channel 
of full width at half-maximum (FWHM) Su, this effecti\^c resolution is approximately 



Selected channel maps from the PPV cube computed with 7„ = —3.5 and 7^, = —3.0 
are shown in Fig. 4. These channel maps are reminiscent of what is observed at 21 cm 
for instance (Joncas et al. 1992). The structures in neighbouring channel maps are clearly 
correlated, with both elongated and diffuse structures. If it is partly due to the channel 
width being subthermal, correlation still exists between more widely separated channels. 



In this section the power spectra of individual channel, integrated emission and centroid 
velocity maps, obtained from the PPV data cubes, are examined and linked to the statistical 
properties of the 3D density and velocity fields. To conduct our analysis, PPV cubes with all 
possible combinations of 3D density and velocity fields were simulated. Fig. 5 and 6 display 
the integrated emission map, the centroid velocity map, 3 selected channel maps - one in the 
line center and two in the line wings - and their associated power spectra, for a subset of 
data cubes. More precisely, Fig. 5 shows the variations of these quantities for 7^, = —4 and 
all values of 7„, while Fig. 6 presents the same analysis for 7^ = —4 and all values of 7^. A 
detailed analysis of these figures follows. 



It has been shown by several authors (see for instance Stutzki et al. (1998) and Goldman 
(2000)) that, for optically thin gas with Gaussian statistics, the power spectral index of the 
integrated emission map is exactly the power spectrum of the 3D density field, under the 
condition d/L> 1 where L is the largest scale on the {x,y) plane and d is the depth in the 
z dimension. A brief summary of the derivation of this result is given in the appendix. This 
is true whether the power spectrum is a power law or not, provided isotropy is maintained. 
This property of optically thin media is powerful as it allows to directly deduce the power 
spectrum of the interstellar 3D density field from the observed map. 
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3.3. 



Channel maps 



4. Statistical properties 



4.1. Integrated emission 
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Our simulations confirm this property for the PPV cube integrated emission (see the 
top panels of Fig. 5 and 6) and for the image of the 3D density field integrated along the z 
axis (see Fig. 7). The fact that the 3D density field and the integrated emission - which is its 
projection onto the plane of the sky - should have the same spectral index does not mean, 
however, that they exhibit the same amount of structure as a function of scale. As already 
suggested in footnote 6, the relevant number in this comparison is the Hurst exponent. Now, 
since 73D = —IHzvi — 3 and ^proj = —'^Hp^oj — 2, the preservation of the spectral index means 
that Hproj = i?3D + 0.5, an effect called projection smoothing. 

On the other hand a cut through the 3D density field would statistically represent the 
original field, and the Hurst exponent would remain the same Hcut = H^j^. The power 
spectrum of such a cut would then have an index ^cut — ~2ifcut — 2 = 730 + 1. An example 
of this is shown in Fig. 7 for a 1-pixel slice of the 3D density cube 

It is also interesting to study the power spectrum of a slice in the intermediate regime, 
where < d < L. To study this, the power spectra of density field shces of increasing 
thickness d were computed. Fig. 8 presents the power spectra of three cases for d—b, 10 and 
20 pixels. The 3D density field used here has 7„ = —3. The power spectra are bent, with 
a spectral index of 7„ at small scales and 7^ + 1 at large scales. The transition occurs at a 
frequency kd — l/2d and provides a method to determine the depth of the observed medium. 
This is in accordance with the pioneer work of von Hoerner (1951) (summarized by O'Dell 
& Castaneda (1987)) who studied the effect of projection smoothing on the index of the 
structure function. Recently, Elmcgrccn ct al. (2001) used this property to determine the 
thickness of the Hi layer in the Large Magellanic Cloud. The absence of such a curvature 
in the power spectra of large scale Hi Galactic emission suggests in turn that the depth 
of the Hi layer sampled by the observations is significantly larger than the map size, i.e. 
much larger than 25 pc in the case of the high latitude cloud observed by Miville-Deschenes 
et al. (2003). This result is also relevant for studies of optically thick media (e.g. molecular 
clouds) where the observed emission comes only from the surface of the cloud. In that case, 
the spectral index measured on the integrated emission map of an optically thick tracer 
should be closer to 7^ + 1. 



^°In this figure, the power spectrum shows a fair amount of dispersion and the index is a bit lower than 
73D + 1. This can be understood by the fact that our sHce through the 3D density cube is not a "perfect" 
one, but indeed has a certain width. The dispersion effect is more pronounced at smaU scales, which are 
comparable to this width. 
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4.2. Centroid velocity 



Another quantity built from PPV data cubes is the centroid velocity map, which, for 
any line of sight, is simply an average velocity weighted by the line temperatures. 



This equation is easily shown to be equivalent to the density-averaged velocity compo- 
nent parallel to the line of sight (Dickman & Kleiner 1985): 



Centroid velocity maps have traditionally been used to study gas kinematics in various 
ISM environments. For instance in the study of the turbulent energy cascade in Hll regions 
(Miville-Deschenes et al. 1995) and in molecular clouds (Miesch & Bally 1994). But one of 
the main difficulties of these analyses is to relate the statistical properties of the centroid 
velocity map to the 3D velocity field. 

To assess this relation in the framework of our fBm simulations, we computed the 
centroid velocity maps from PPV cubes for every combination of 7„ and 7^,. Some of them 
are shown in Fig. 5 and Fig. 6, along with their power spectra, and all measured spectral 
indexes are given in Table 1. In the range of density and velocity indexes tested, the spectral 
index 7^7 of the centroid velocity map is remarkably equal to 7^, whatever the power spectrum 
of the density field. Fig. 5 illustrates this behaviour: the velocity centroid maps computed 
for the different values of 7^ are the same, and their power spectra all merge into the same 
spectrum. 

This result is easy to obtain analytically for a uniform density, as is shown in the 
Appendix (see also Dickman & Kleiner (1985)). But, to our knowledge, this is the first 
time that it is demonstrated for stochastic density fields. Fig. 9 displays the full scatter 
plots of power versus wavenumber for several (7^, 7^) combinations. It appears that the 
dispersion around the power law decreases with steepening density spectra, and in the limit 
of a uniform density field, isotropy is recovered. Conversely, for a given 7„, the centroid 
spectra are more scattered for steeper values of 7^. This may be explained by the fact that 
density fluctuations manifest themselves more strongly with smoother velocity fields. 
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T^u^uix, y,u)5u 
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Further tests of this result have been performed. For instance, in accordance with 
Eq. 10, it does not depend on the gas kinetic temperature, with 7c being equal to 7^ in all 
cases considered. More interesting was the test of non-fBm density and velocity fields. One 
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of these involved a velocity field for which the power spectrum was bent, varying from a 
behaviour at large scales (small k) to a behaviour at small scales (large k). For this, two 
fBms with respective indexes -3.0 and -5.0 were built and their Fourier amplitudes added with 
proper weights^^. After introduction of a random phase map, an inverse Fourier transform 
was performed. The velocity field computed this way was then used with two density cubes, 
one with 7„ = —4.0 and the other being the uniform density cube, to simulate two PPV 
data sets. Velocity centroids for both cases were then computed and their power spectra are 
shown in Fig. 10. For the uniform density case, unsurprisingly, the power spectrum of the 
velocity centroid has zero dispersion and exhibits exactly the same variation as that of the 
full 3D velocity field. For the 7^ = —4.0 case, despite the dispersion, the behaviour remains, 
showing the consistency of the method, even for isotropic velocity fields with non-power-law 
power spectra. 

As a second extension of the types of fields tested, density cubes with increased fluctu- 
ations were built by exponentiation of fBm fields, resulting in lognormal distributions. The 
original fBm field used had an index 7 = —4.0 and the resulting density values ranged from 
~ 1 cm~^ to ~ 2500 cm~^, with a mean of ~ 90 cm^^. This density cube was used in con- 
junction with all of our fBm velocity fields to compute PPV data sets and velocity centroids, 
the power spectra of which are shown in Fig. 11. For such large density fluctuations, the 
dispersion is noticeably increased with respect to the fBm case, but the main result that 
7c = 7?, holds when averaging in wavevector space. 

Considering the possibility that the complete agreement between 7^ and 7,u, without 7„ 
entering into the picture, could be caused by the absence of density- velocity correlations, we 
chose to compute 3D density and velocity fBm fields featuring a certain level of correlation. 
The way we did this is based on the idea that most, if not all, of the spatial structure of 
images is contained in the structure of their Fourier phases. Starting from this, we introduced 
correlations between two fBms by setting a number of phase points to the same value in both 
fields. The points were chosen randomly, and their number varied between 10% and 90% of 
the total number of pixels. The amplitudes were computed separately. 

Once again, we found that the velocity centroid spectral index was exactly the same 
as that of the three-dimensional velocity field, whatever the density spectral index, and 
whatever the amount of correlations. From this remarkable result we infer that velocity 
centroids are robust observables for recovering 3D velocity statistics from optically thin 
lines. This result has obvious applications for studies of ISM statistical properties and of 



^^This is to ensure that the velocity field is not dominated by either behaviour over the whole range of k, 
in which case the curvature of the power spectrum would be imossible to see. 
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interstellar turbulence (see Miville-Deschenes et al. (2003)). 

We are aware that the method used to introduce density-velocity correlation in our 
simulations does not comply with (magneto-)hydrodynamics, as for instance high densities 
would correspond to high velocities along the line of sight, but this technique allows us to 
test easily various correlation degrees while staying in the fBm scheme. Obviously, it would 

be useful to make a similar analysis on physically relevant fields, such as MHD simulations, 
to study the effects on the centroid velocity power spectrum of physical density-velocity 
correlation, sharper density/velocity contrasts and non-Gaussian statistics, and thus assess 
the robustness of this result in more realistic conditions. 

One should also keep in mind that our results may not be generally applicable to all 
interstellar conditions. For instance the statistical properties of centroid velocity maps can 
be significantly affected by systematic motions like rotation or expansion, if proven not to 
be part of the turbulent cascade itself. In this case, such systematic velocity fields should be 
filtered prior to any analysis of turbulent motions (Kleiner & Dickman 1985). 

4.3. Channel maps and velocity slices 

In a recent study, LPOO showed that, for Gaussian random fields, the spectral index of 
the 3D velocity field can be determined by studying the variation of the spectral index of 
velocity slices^^ as a function of their thickness AV. Subsequently, Lazarian et al. (2001) 
and Esquivel et al. (2003) analyzed and confirmed this proposition using MHD simulations. 
In both these works the power spectrum of the simulated velocity field is forced to be a 
power law (a method called spectral modification) to extend the limited power law range 
of MHD simulations. This method preserves the phase information and thus much of the 
expensively-obtained physics. We propose here to also test the theoretical predictions of 
LPOO on fBm simulations, because they have intrinsically well defined power spectra and 
obey Gaussian statistics, the working hypotheses of LPOO. 

4-3.1. Thickness of velocity slices 

The main prediction of LPOO deals with what they call thin and thick velocity shoes of 
the PPV data cube. These properties are scale-dependent: A velocity slice of width AV is 
considered thin at a scale I if the velocity dispersion in the observed medium at that scale 



'A velocity slice of thickness AV is the integrated emission of the PPV cube between v and v + AV. 



a{l) is larger than AV. Conversely, a velocity slice is said to be thick at scale I if a{l) < AV. 
Since the velocity dispersion depends on I according to 




(-7^-3)/2. 
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where (t{L) is the total velocity dispersion^'^ (set to 3 km s~^ in our simulations), velocity 
slices arc thinner at large scales than at small scales (7^ < —3). In fact, from Eq. 11, the 
smallest scale Ithin for which a velocity slice of width AV can be considered thin is given by: 



Given the values of AV and cr(L) sampled in our simulations, the velocity slices under study 
can be thin or thick at all scales, in which cases they are respectively said to be fully thin 
or fully thick. They can otherwise exhibit a transition between the two regimes, being thin 
at large scales and thick at small scales. Obviously, a velocity shoe will be fully thin if all 
sampled scales are larger than Ifhin, and fully thick if the largest scale in the map is smaller 



At that point it is important to check whether or not our simulations allow us to 
investigate the thin regime. The velocity dispersion at the 1 pixel scale (see Table 2) is given 

by Eq. 11, and is to be compared with the smallest AV available, which is determined by 
the effective spectral resolution (sec Eq. 8), namely Sve/f = 1.4 km s^-*^ for 6u = 0.25 km s~^ 
and T = 100 K. It appears that, for 7^,, = —3, there are meaningful velocity slices (with 
AV > Sveff) in the fully thin regime, since a{l =1 pixel) > Sv^ff. For smaller values of 
7„, however, any velocity slice will be thick at small scales. Nonetheless, the values of khin, 
given for all 7„ in Table 2, show that a AV — Svgff velocity slice is thin on a large range of 
scales, even for 7^ = —4.5. In order to have access to fully thin velocity shces for all 7^, an 
increased dispersion a{L) or a lower temperature would have been necessary (see Eq. 11), 
but that would have led to more velocity channels in the PPV cube^^ and therefore to more 
shot-noise, which has a dramatic effect on simulated channel maps, as will be discussed later. 

The condition expressed above for a velocity shoe of width AV to be fully thin calls 
for a remark. The criterion a{l) > AV is given for two-dimensional scales I in the plane of 
the sky, but the velocity dispersion of the gas sampled at that scale by the observations is 
actually a function of the 3D scales ^30. Many of those, when projected, give the same 2D 



2 




(12) 



than I 



'thin- 



L is the largest scale in the image. 

The spectral resolution 6u being kept constant, in order to keep Svgff unchanged. 
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scale I, but their velocity dispersions can be quite different, and an average should enter in 
the thinness condition. Its present form however ensures that, since Z3D > I, if the velocity 
slice is thin at scale I, it will also be thin at all unprojected scales ^30, according to Eq. 11. 

According to LPOO, fully thin velocity slices, when integrated, produce maps whose 
spectral index is (for 7„ < —3) 

Ithin = -3 - ^''^ , (13) 

and fully thick velocity slices exhibit a spectral index 

^thick = -3 H — . (14) 

When a substantial part of the emission is integrated over, the slices are said to be very thick 
and their spectral index tends to that of the 3D density field, as shown in 4.1. 

To test this prediction, the spectral indexes of velocity shces of increasing width were 
computed, for all our PPV data cubes, and at a temperature of T —100 K. The power 
spectra of selected velocity slices with difi^erent widths AV are shown in Fig. 12. Each 
plot corresponds to a different set of (7^, 7^,) and shows, for the given PPV cube, the power 
spectra of four velocity slices (open circles) with AV = 1.25, 1.75, 2.25 and 10.0 km s"^(from 
bottom to top). In each case, the solid lines represent the predictions of LPOO for the thin 
and thick regimes (see Equations 13 and 14), and the small vertical marker indicates the 
transition scale Imn, given by Eq. 12. The velocity slice under consideration is thin for 
scales on the left of this mark, and thick for those on the right. The power spectra shown in 
Fig. 12 are in general agreement with the predictions of LPOO, given the selected values of 
AV^, 7„ and 7^. Bent power spectra are observed, with a break close to Imn- However, the 
narrow range of scales over which the velocity slices are thin makes the comparison with the 
predictions uneasy. Furthermore, the power spectra of velocity slices for 7„ = —4.5 move 
away from the LPOO prediction at small scale, because these velocity slices are becoming 
very thick at those scales. This behavior is indeed correlated with 7^. Another feature of 
Fig. 12 is the increase of power at small scales for 7^ = —3, due to shot-noise. 



4.3.2. Finite-size or "shot-noise" effect 

The computation of PPV cubes from finite size 3D velocity fields leads to an artificial 
increase of high-frequency fiuctuations in channel maps (Lazarian et al. 2001; Brunt & Heyer 
2002a; Esquivel et al. 2003). As pointed out by Esquivel et al. (2003), the limited size of 
the simulated fines of sight introduces sharp edges and almost empty channels in the PPV 
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cube. This effect is more important witfi a small effective spectral resolution. Such a high- 
frequency flattening is seen in the power spectra of channel maps for 7^, = —3 in figures 6 
and 12. As expected, the amplitude of this high-frequency component increases when the 
size of the 3D fields decreases. This effect biases the determination of simulated channel 
maps power spectra. 

One sohition to reduce the effects of shot-noise is to apply a Gaussian smoothing to 
the PPV spectra, simulating the effect of thermal broadening. But such a smoothing also 
increases the effective spectral resolution (see Eq. 8) which results in a reduction of the 
accessible thin regime. An example of this is shown in Fig. 13 where we plot the power 
spectra of velocity shces with AV = Sveff, for two PPV cubes ((7^ = —4, ■jy = —3) and 
(7„ = —3, 7^ = —4.5)) computed with different temperatures {T — 1, 10, 100, 200 and 500 
K). The 7„ = —3 case is particularly interesting as velocity shces are thin at all scales as 
long as AV < a{L) (see Eq. 11), which is the case here {(j{L) — 3 km s~^ and 5veff — 2.9 
km s~^ for the highest temperature T = 500 K). This allows for a separate study of the 
effect of Gaussian smoothing on shot-noise reduction and velocity slice thinness. The top 
panel of Fig. 13 shows that Gaussian smoothing significantly dampens shot-noise. Up to a 
temperature of 500 K the power spectrum steepens through reduction of the shot-noise effect, 
but it seems that, even at this temperature, shot-noise has not been completely removed. 
At higher temperatures, the condition Sveff < o-{L) would not be met any longer. 

The bottom panel of Fig. 13 illustrates the effect of Gaussian smoothing on both the 
shot-noise and the thickness of the velocity slice. At low temperature (T =1 K), the velocity 
slice is thin on almost all scales but the effect of shot-noise can be observed at high-frequency 
as the power spectrum is above the prediction of LPOO. As expected, the shot-noise amplitude 
and the range of scales at which the velocity slice is thin both decrease with increasing T. 
At high temperature, the shot-noise is completely removed but the velocity slice is fully 
thick. This shows how the use of Gaussian smoothing is necessarily a compromise between 
the reduction of shot-noise and the thinness of velocity slices. From the power spectra 
shown in Fig. 12, shot-noise appears to be significant only for 7^, = —3, indicating that the 
simulation parameters selected (cr(L) = 3 km s~^, T = 100 K and 6u = 0.25 km s^^) provide 
a good compromise, for most of our 7^ values, between shot-noise reduction and channel 
map thinness. 

4-3.3. Velocity slice analysis 

The method proposed by LPOO to determine 7„ and 7^, is to look at the spectral index 
ilsiice) of velocity slices of increasing width A^. According to LPOO, the curve of ■jsHce 
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vs should reach two asymptotic values at small (thin) and large (very thick) AT/ that 
allow the determination of 7„ and 7^. This velocity slice analysis has been performed on 
our simulations, and the results are presented in Fig. 14, which shows the variations of the 
index '-jsUce (fitted over the whole range of k) as a function of the velocity width AV. The 
curves, one for every set of (7„, 7^,), are grouped by their 7^, value. The shaded areas in 
this figure indicate the fully thin, fully thick and very thick regimes but also the transition 
regime where the velocity shce is thin at large scales and thick at small scales. The vertical 
dashed line indicates the effective spectral resolution of our simulations. The plateau seen 
on the left of this fine reflects the fact that velocity structures cannot be resolved due to the 
flnite spectral resolution. 

As Ay enters the very thick regime, the spectral index ^sUce converges to the 7„ value, 
in agreement with the fact that the integrated emission has the same power spectrum as 
the 3D density held. At the other end, we confirm the prediction of LPOO that the power 
spectrum of narrow velocity slices is independent of the density field. For a given 7^, all curves 
converge to the same value for small AK, except in the case of 7^,, = —4.5, as explained in 
the next paragraph. This strong correlation between the 3D velocity field and channel maps 
is also seen in Fig. 5 and 6. First, in Fig. 5, the channel maps' structure is dominated by the 
velocity field as it varies very little with 7^. Similarly, the channel maps power spectra for a 
given 7^ are the same for all values of 7„. The influence of the velocity held on the channel 
maps structure is also manifest in Fig. 6 where channel maps power spectra vary with 7^. 

Back to Fig. 14, the curves for 7^ = —3 rise at small AV, up to a value of -1.8. As the 
spectral index '-jsUce has been computed over the whole k range, this trend can be completely 
attributed to shot- noise, which has a significant effect for 7„ = —3 (see Fig. 13). For 
7^ = —3.5 the effect of shot-noise is reduced, and in that case, the curves are in relatively good 
agreement with the predictions of LPOO. There is a significant plateau in the thick regime, 
at a value very close to the predicted ^tuck value (see Eq. 14). Moreover, for Ay = ^veffi we 
are very close to the fully thin regime and the spectral index measured there is consistent 
with LPOO. For 7^, = —4 and -4.5, the concordance with LPOO is not that clear. In these 
two cases the effect of shot-noise is completely removed by the Gaussian smoothing but, on 
the other hand, we are very far from the fully thin velocity slices. Some kind of plateau is 
observed in both cases within the thick regime, but not at the predicted thick value. This is 
due to the fact that the smallest scales are actually becoming very thick. This is why, in the 
case of 7^ = —4.5, the curves are dependent on 7„ and do not converge to the same value. 
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4.3.4- Discussion 



Except for 7t, = —3.0, we have seen that shot-noise has a hmited effect on our simulations 
because of the Gaussian smoothing apphed. But this is at the expense of a narrowing of 
the accessible thin regime: fully thin velocity slices with AV" > Svgff are present only when 
— —3.0, In these conditions, the interpretation of the ^sHce vs AV curves (see figures 13 
and 14) is ambiguous, which makes the determination of 7^ very difficult. The power spectra 
of fully thick velocity slices can actually depend on 7^, which makes it impossible to deduce 
7^ from ^thick using Eq. 14. Therefore, due to the relative thickness of velocity slices in our 
simulations, we could not demonstrate whether the LPOO method can be used to determine 



One practical hmitation of the LPOO method, however, is that it requires to find out 
whether or not the fully thin regime is reached, for a given observation. If a significant 
asymptotic plateau is found at small AV in the 'jsiice vs AV curves, one might think that 
it is the case. But unfortunately, as seen in Fig. 14, such a plateau could also be observed 
in the thick regime, and this could lead to a wrong estimate of 7^,. In addition, the only 
way to estimate the position of the transition scale khin between the two regimes for a given 
velocity slice is by using Eq. 12, which requires 7^. Therefore, there is no independent way 
to determine whether thin velocity slices have been reached or not. In that respect, the 
determination of 7^ in the Small Magellanic Cloud by Stanimirovic & Lazarian (2001) might 
be questionable as they do not seem to reach the thin asymptotic regime. 

It is important to estimate in what conditions LPOO's method could be used on real 
observations, for a typical high latitude Hi cloud for instance. We have shown that the nar- 
rowest relevant velocity shce for this analysis has a width of AV — Syeff ~ •\/2.16 /csT/m 
(here we consider that the effective spectral resolution is dominated by the thermal broad- 
ening). But to estimate 7^ using LPOO's method one must have a plateau at small AV 
in the (JsUce vs AV curve, which calls for fully thin velocity slices on a significant range of 
AV, lets say for 6yeff < AV < 2 6veff- We recall that velocity slices are thin at scales 
i > hhin where cr{lthin) = AV. Therefore to evaluate on what scales LPOO's method could 
be used to deduce 7^,, one must assess what is Ithin for a velocity shce of AV — 2Sveff- As 
shown by Miville-Deschenes et al. (2003), the energy spectrum of Hi is compatible with the 
Kolmogorov (1941) cascade {a{l) — Af -^^). In these conditions, khin for a velocity shce of 
AV^25yeff is 



The normalisation factor A can be estimated from high-latitude Hi observations like in Jon- 
cas et al. (1992), where the internal velocity dispersion at a scale of 1 pc is of the order of 4 




(15) 
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km s~^, significantly larger than that observed in molecular clouds (Larson 1981). If we con- 
sider Cold Neutral Medium (CNM) gas at T=100 K, the smallest thin scale Ithin for a cloud 
at a distance of 100 pc is ~ 8', which is significantly larger than the angular resolution (1') of 
the 21 cm radio-interferomctcr at Dominion Radio Astrophysical Observatory (DRAO). Fur- 
thermore, if we consider the presence of warmer gas (Warm Neutral Medium (T = 5000-8000 
K) or thermally unstable Hi (T ~ 1000 K)), the thermal broadening would increase Ithin 
(which scales as T^/^) and make the determination of 7^ using the LPOO method only possible 
at scales larger than a few degrees, at least for Hi in the solar neighborhood. In this respect, 
the velocity centroid method is better suited to determine the 3D velocity power spectrum 
as it does not depend at all on the gas temperature and on the thickness of the velocity shces 
in the observation. Therefore the velocity centroid method allows the determination of the 
velocity power spectrum down to the angular resolution of the observations. 



5. Conclusion 

In this paper three-dimensional fBms with various spectral indexes were used to simulate 
interstellar 3D density and velocity fields. Prom these we computed several spectro-imagery 
observations assuming optically thin line emission. Our results confirm that the spectral 
index 7^ of the integrated emission is the same as that of the 3D density field 7„, when the 
depth d of the medium probed along the line of sight is at least of the order of the largest 
transverse scale observed. For smaller depths, the power spectrum of the integrated emission 
shows two asymptotic slopes of 7„ and 7„ -|- 1 at small and large scales respectively, with 
the transition occurring at a frequency k — l/2d. This property of the power spectrum of 
integrated emission could be used to estimate the depth of interstellar clouds and to infer 
the 3D statistical properties of optically thick media. 

The main result of our study is that the spectral index of the maps of hne centroid 
velocities is the same as that of the 3D velocity field. This result is particularly robust because 
it holds for any gas temperature and whether density- velocity correlations are present or not. 
This similarity is also recovered if the power spectrum of the 3D velocity field is not a power 
law or if the density fluctuations follow a lognormal distribution. This remarkable property 
of optically thin media may be used to estimate the 3D energy spectrum of interstellar 
turbulence in various ISM environments. 

We also compare the power spectra of velocity slices with the predictions of LPOO. Due 
to the limited resolution of our simulations we arc not able to confirm their predictions. On 
the other hand we show that the method proposed by LPOO is in practice barely applicable to 
real observations, and we suggest that a more robust determination of 3D velocity statistics 
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lies with velocity centroids. 



A. Spectral index of the integrated emission 

In our case where the medium under study is optically thin, the integrated emission 
map is simply the sum of all velocity channels: 

W{x, y) = Y, ^u{x, y, u)Su. (Al) 

u 

Using Eq. 6, reverting the order of summations over z and u, and with our values of 6z and 
Su, this translates into 

W{x, y)^aY^ n{x, y, z)5z (A2) 

2 

where a is a numerical coefficient coming from the unit conversion. The important point is 
that W is proportional to the sum over z of the 3D density cube, which by construction is a 
fBm with spectral index 7„. To compute the spectral index 7vi/ of the integrated emission, 
we consider the limit for which the summation can be replaced by an integral over z 

W{x,y)—a J n{x,y,z)dz. (A3) 

Writing the Fourier transform of the 3D density as a function of the three-dimensional wave 
vector k, we have 

n(k) = /// n{r)e-^'''^-''dxdydz. (A4) 



Which becomes, if we consider the kz — cut through fi(k) and integrate first over z, 

n{k^, ky, = JJ ^^^'^^ e-^^"(^-^+^^^Ma:dy (A5) 

meaning that the Fourier transform of W is proportional to the kz — cut through h. We 
therefore have 

\W{k,,ky)\' (X \h{k,,kyM' OC {kl + kly"/', (A6) 

and it is then easy to see that Pw and P„ follow the same law, and thus 

'fw = In- (A7) 

This calculation can be fully duplicated to compute the spectral index of the centroid velocity 
map, in the case of a uniform density field. In this situation, the centroid velocity map 
(Eq. 10) is just 

C(x,y)^(3 J v{x,y,z)dz. (A8) 
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We therefore have a similar result which reads 



7c = Iv- 



(A9) 
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Fig. 1. — Typical two-dimensional 257 x 257 fractional Brownian motion image with its 
power spectrum, which by construction follows a power law with an index equal to -3.5. 
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Fig. 2. — Typical fractal line of sight extracted from the 3D density (bottom, 7„=-4.5) and 
velocity (top, 7„=-3.0) fields. 
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Fig. 3. — Synthetic spectra for the hne of sight described in Fig. 2, at kinetic temperatures 
T = 1 K, 10 K, 100 K and 500 K. The spectral resolution was set to Su = 0.25 km s^^ 
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Fig. 4. — Selected channel maps from a typical simulated PPV data cube, shown here for 
comparison with actual observations. In this case the power spectrum index of the density 
and velocity fields are respectively 7„ = —3.5 and 7^ = —3.0. The gas temperature used is 
T — 100 K (thermal velocity dispersion ath — 0.91 km s~^) and the spectral resolution is 
5u — 0.25 km s~^. The channel number is indicated in the upper left corner of every channel 
map. 
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Fig. 5. — Integrated emission, centroid velocity and three velocity channels for 7^, = —4 
simulations. The four image columns correspond to four values of 7„ specified above each 
column. The power spectra of all 2D maps shown here are displayed on the far right column, 
the solid, dashed, dotted and dash-dotted lines representing respectively 7„ = —3.0, -3.5, -4.0, 
and -4.5. This figure highfights the fact that the power spectrum of the integrated emission 
is completely determined by the 3D density field (we verify that their power spectrum index 
are equal) . But the most striking feature of this figure is the fact that the centroid velocity 
is independent of the density field and that the power spectrum of the centroid velocity is 
equal to the power spectrum of 3D velocity field. In addition, the structure of channel maps 
and their power spectra are mostly determined by the velocity field (we observe a slight 
variation of the channel maps structure with different density fields) . 
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Fig. 6. — Integrated emission, centroid velocity and three velocity channels for 7„ = —4 
simulations. The four image columns correspond to four values of 7^, specified above each 
column. The power spectra of all 2D maps shown here are displayed on the far right column, 
the solid, dashed, dotted and dash-dotted hnes representing respectively 7^, = —3.0, -3.5, -4.0, 
and -4.5. This figure confirms the fact that the integrated emission is completely independent 
of the velocity field (top row) and that the power spectra of the velocity centroid and channel 
maps are completely dominated by the 3D velocity field. This figure also shows the effect of 
"shot-noise" discussed in the text that produces a flattening at small scales of the channel 
maps power spectrum, for shallow (7^ = —3) velocity fields. 
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Fig. 7. — Map and power spectra of the fully integrated emission (top row) and of a single 
slice (bottom row) through one of the 3D fBm density cubes, in this case 7^ = —3.5 and 
Su^ 0.25 km s-^ 
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Fig. 8. — Power spectra of integrated emission slices of a 7„ = —3.0 density cube. The 
depth of the shces are ci = 5, 10 and 20 pixels, from left to right. The solid line is the power 
spectrum of the integrated density shoe. The dotted lines indicate —3.0 and —2.0 slopes and 
the vertical dashed line shows the l/2d scale. 
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Fig. 9. — Power spectra of centroid velocity maps for various (7„, 7^) combinations. Each 
subplot corresponds to a given 7^, and features the cases for 7„ = —3.0 (black scatter plot), 
7„ = —4.0 (grey scatter plot), and uniform density (thin black line). The latter is actually 
a scatter plot with virtually zero dispersion, as expected. 
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Fig. 10. — Power spectra of centroid velocity maps constructed from a velocity field with 
a curved power spectrum (see text). Plots for the uniform density case (black) and the 
7„ = —4 case are presented, along with the expected asymptotic behaviours at small {k~^) 
and large (k^^) scales. 
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Fig. 11. — Power spectra of ccntroid velocity maps constructed from a lognormally dis- 
tributed density field (see text). Each subplot corresponds to a given 7^,, and features the 
lognormal density case (black scatter plot), and the fBm density case with 7^ = — 4(grey 
scatter plot). The black line represents the expected behaviour in k~"'^ . 




Fig. 12. — Power spectra of velocity slices for all (7^, 7„) combinations. On each plot, the 
open circles represent the power spectra of velocity slices for AV = 1.25, 1.75, 2.25 and 
10.0 km s^^ (from bottom to top). The solid line on each curve is the prediction of LPOO 
(see equations 13 and 14), normalized to fit around the thin-thick transition (indicated by 
a vertical marker). On each curves, the scales are thin (thick) on the left (right) side of this 
vertical line. 
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Fig. 13. — Power spectra of 1 pixel velocity slices of width AV = 1.25 km s~^ for (7^ = —4, 
7„ = —3) (top) and (7^ = —3, 7^, = —4.5) (bottom). On each plot, the five curves (open 
circles) represent PPV cubes computed with T = 1, 10, 100, 200 and 500 K. The solid line 
on each curve is the prediction of LPOO (see equations 13 and 14), normalized to fit around 
the thin-thick transition (indicated by a vertical marker). The slices are thin for scales on 
the left side of these vertical fines, and thick for scales on their right side. 
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Fig. 14. — Velocity slice spectral indexes as a function of velocity width AV, computed for 
all (7„, 7^) combinations. Curves are grouped by their 7^, value. The white and grey zones 
indicate the "thin", "transition", "thick" and "very thick" regimes. The thin regime is for 
AV^ < (t{1 = 1 pixel). The thick regime is for AV > cr{L) (the velocity width is larger than 
the velocity dispersion at the largest scale). The transition regime corresponds to velocity 
slices with intermediate AV^ (they are thin at large scales and thick at small scales). The 
very thick regime is where '^siice = In {AV ~ 20 km s^^). The vertical dashed line indicates 
the effective spectral resolution of our simulations (()t,e// = 1-4 km s~^). The horizontal 
dotted lines indicate LPOO's predictions for ^thin and ^mck (see equations 13 and 14). 
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Table 1: Spectral index of the centroid velocity {jc) for the simulated PPV data cubes 
computed in the analysis. The uncertainty of the measured spectral indexes is 0.1. 



7« 


(t(/=1 pixel) 


khini^V = 6Veff) 


-3.0 


3 km 


pixel 


-3.5 


0.9 km s"^ 


4 pixels 


-4.0 


0.3 km s-^ 


24 pixels 


-4.5 


0.08 km s-^ 


42 pixels 



Table 2: These values were computed for T = 100 K, Su = 0.25 km s ^, (j{L) = 3 km s ^ and 
L=129 pixels. The effective spectral resolution of our simulation (see Eq. 8) is 1.4 km s~^. 



